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Abstract. We consider the spin-glass phase of the Sherrington-Kirkpatrick model in the presence of a 
magnetic field. The series expansion of the Parisi function q{x) is computed at high orders in powers of 
T = Tc — T and H . We find that none of the Parisi-Toulouse scaling hypotheses on the q{x) behavior 
strictly holds, although some of them are violated only at high orders. The series is resummed yielding 
results in the whole spin-glass phase which are compared with those from a numerical evaluation of the 
q{x). At the high order considered, the transition turns out to be third order on the Almeida-Thouless 
line, a result which is confirmed rigorously computing the expansion of the solution near the line at finite 
T. The transition becomes smoother for infinitesimally small field while it is third order at strictly zero 
field. 

PACS. 75.10.Nr Spin-glass and other random models - 02.30.Mv Approximations and expansions 



Mean-field theory is commonly used as a starting point 
(i.e. zeroth order approximation) in the perturbative treat- 
ment of finite dimensional, short-ranged systems. As such, 
it is usually very simple, almost trivial, in ordinary sys- 
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terns like the ferromagnet. On the contrary, the glassy 
phase of disordered systems is highly non-trivial even in 
mean-field theory: in the Ising spin glass, the prototype 
of such systems, the thermodynamics of the low tempera- 
ture phase is built on an ultrametrically structured order 
parameter matrix proposed by Parisi (see [Ij and some of 
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the reprints from this book). Although the region close to these are the only ones which may be exact. The order 
the zero magnetic field critical point is easily accessible, of the transition when crossing the AT-linc can also be 
to get quantitative results for physical observables deep deduced from the projection hypothesis, and it turns out 
inside the glassy phase is a rather difficult task. to be second order in the Ehrcnfest sense, at least in finite 

field (The magnetic susceptibility has, for example, a 

The survival of the mean-field picture, at least quali- 

finite, albeit small jump.) Nevertheless, subsequent calcu- 

tatively, in finite dimensional, short-ranged systems has 

lations and arguments [7lllU| seemed to provide evidence 

been a controversial and long debated problem in the 

that the transition is third order, just as in zero magnetic 

last decades. It is clear that a deeper, and more precise, 

field. Recent numerical simulations on the SK model in a 

understanding of the mean-field thermodynamics of the 

field are discussed in 

spin glass phase is a prerequisite to make any compari- 

. 1 1 i- 1 -ii 1 • 1 i rni In this work we study the Parisi solution of the Sherrington- 

son, numerical or analytical, with physical systems, ihe ° 

Ti • • m 1 /n rn\ 1 ,1 r r, £ Tji r ^ 1 Kirkpatrick (SK) model in the presence of an external field 

Parisi- loulouse (Pal) hypotheses of Ret. |4| (also known i v / f 

. 1 ,1 N _L 1 iz i i • through the methods developed in 1121 . In particular, the 

as projection hypotheses) were not only a first step m ° 

^1 . ,. ^. 1 i 1 1 , 1 1 , r • 1- -i computation of the high order expansion of (7(a;) in power 
this direction, but had the advantage of simplicity, a com- ^ ^ ^ i\ / t- 

. s- £!ii_Li A i of the reduced temperature and of the field allows us to 

mon feature of mean-held theory. An extension to systems 

.^1 s- i- • i r n 1 finally answer the problem of the validity of the PaT hy- 

with an average ferromagnetic interaction followed soon j j j 

nn 1 1- - ■ r ;i J potheses 011^1121 ■ The PaT hypotheses are a set of scaling 

|5j, whereas a scaling assumption tor the order parame- ^ jf a 

r i- I \ i r ^ ■ rn Ti relations concerning the behavior of the function qix) un- 

ter function q[x) was put forward m p|, see eqs.|iiJ It was ' 

n . -ij 11 i-n ii der changes of the external fields. As we shall briefly recall 

first pointed out m ret. 6 that the whole set of hypotheses ° 

i-iii^- -i A i below, they allow to compute all the quantities of interest 

cannot be exact, as it leads to inconsistency. An argument ; j i i 

, , . r iwi ii i ii ■ i- 1 ii • of the spin glass phase without solving the Parisi equa- 

was presented m ret. |7| that the projection hypothesis t- a t- o -± 

r ,1 i- i- i J- i j_i • • 1 tions but simply using the SK replica symmetric solution 

tor the magnetization contradicts the maximization rule i j a j- j 

r ^, r ,1 _L ii-ii-i-i Ai -1 evaluated at the AT line 8 . It turned out from our analy- 

tor the free energy, at least near the instability, Almeida- - - 

mi 1 /Arn\ 1- m A 1 • Til i Tir • sis that no one of the PaT statements strictly holds, some 

ihouless (Ai), line |8|. An analysis of the exact, dittusion- 

T^ j.rr ,-1 ,• J ■!• ,1 ,1 1 • of them, however, are so highly consistent with the numer- 

like dinerential equations describing the thermodynamics ° 

r , 1 £^ ij 1 1 i J • A TT ical estimates of the q(x) that their approximate nature 

ot the mean-field glassy phase was presented m j9i. Us- ' 

1 1 1 _L • -i can be detected only at high order in the series expansion, 

mg series expansions, with relatively short series, it was jo 

concluded in this paper that the order parameter function If one evaluates numerically the q{x), one finds that 

q{x) and the breakpoint xi depend only on temperature, the magnetization seems to be independent of the tem- 

and, among the set of statements of the PaT hypotheses, perature in the spin glass phase; similarly the self-overlap 
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qea seems to be independent of the magnetic field. This laws: 
motivated the PaT hypotheses which assume this fact to 
be exactly true, i.e. 

M{H,T) = M{H), qEA{T,H) = qEA{T). (1) 



qoiH,T) ^qoiH)= qsKiH,TATiH)), 
q[x,T,H) =q{x,T)^qu{x/T), (6) 
qEA{T,H) = qEA{T). 



The universal temperature independent function q^ {y) can 

As a consequence, the free energy is additive, i.e. we have , ; j <• ;i i .• ri ^ ^ j i i • i 

: l^g computed from the relation q{x)dx = 1 — j , which 

F(T,H) = Fi{T) + F2(H). Furthermore, the energy is , , , . 4.- -c u j ^i, i i j t frr\ 

^ ' ' L\ / ^\ / I oj holds m zero magnetic held, and the knowledge of g£;yi(J ). 

also additive, while the entropy depends only on the tern- r^, „ r- . . i i ; • i ■ 

i heretore any quantity of interest can be determined sim- 

perature. The previous equations show that if we want , ... ^ Arn t i ot^ i ^■ 

ply projecting it from the Ai line, where the bK solution 

to compute M(H), we can compute the magnetization . ^ , • ^ i- ^ i i /t^ot-)\ 

^ ^' ^ ^ IS correct, down into the replica symmetry broken (RSfij 

at any temperature in the spin glass phase provided the , r • n rn i i i 

^ o jr- jr- phase; for this reason the Pai hypotheses are also known 

field is equal to H\ in particular, it can be computed at ,, . , „ t ^ •, ^ i ■• i 

as projection hypotheses . fn it was already noticed 

the boundary of the spin glass phase, i.e. on the AT line, ,i ; ;i n rp i ,1 u iu- 

i o ; I th&t the Pai hypotheses must be wrong somehow, this 

where the replica symmetric solution is the correct one: i 1 1 • i- • ^1 1 r i-n- 

was deduced indirectly comparing the values of different 

expressions for the specific heat. Instead, here we are able 

M{H)^M{H,T)^Msk{H,Tat{H)). (2) 

to check the relations (j^Jl directly. 
Here the function Hat{T), or its inverse Tat{H), param- The equations determining q{x) follows from Parisi dif- 



eterizes the AT line. The argument can be generalized to ferential equation 2 and the stationarity of the free en- 
any quantity that depends on either the temperature, or ergy. They have appeared many times in the literature, see 
the magnetic field alone. For instance, we have ref. [HltTTIIMn] . therefore we choose don't report them 

here. Besides the function q{x) they involve two auxiliary 



SiT) = Sir, H) = Ssk{T, Hat(T)), (3) 
qEA{T) = qEA{T,H) ^ qSK{T,HAT{T)). (4) 



functions P{x,y) and m{x,y). Introduced as mathemati- 
cal tools to perform the computation, these functions have 
also a physical interpretation: the function 'm{x, y) repre- 
Furthermore, if we consider the expression for the energy ^^^^^ magnetization of a given site properly averaged 



\ / \ (see chap. V in Ref. [li) on a cluster of states with mutual 

E^-—[l~jq\x)dx\-MH, (5) 

^ overlap q{x) in presence of a frozen effective field y, accord- 



we see that a possible way to have an additive charac- ingly m(l,?/) — tanh(/3?/); the function P{x,y) represents 
ter is to guess proper scaling laws also for qQ{H,T) and the probability distribution function over the disorder of 
q{x,H,T). This is at the origin of the full PaT scaling the effective field y, in particular P{l,y) represents the 
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distribution of the cavity field 1 . These functions have 
also received a dynamical interpretation in |13| . 

We have solved the equations by series expansions around ^ 22 256r 3808r^ 



q{x,T,H) ~ q{x,T,0) 
P 



H = and T = Tc = 1. In the following we report the 



„6 



power series of the quantities of interest at lowest orders. ^ \ 9 3 / 35 ^ 

To make easier the comparison between order of magni- 



5 15 45 
-64 64r\ 2144 , 

132 88/ 448 p^r 



5 3 

tude we write down the series expansion in terms of the 312/ a;^ r^/o\ 

H hC(8j, (yj 

variables r = 1 - T and p = (3/4)i/3i?2/3 ^hey are in- 

1 r t2 t3 

tended to be of the same order of magnitude such that the f(T, H) = — (1 — T)ln2 \ 

^ ■'^ ' ^ ^ ^ 4 4 4 12 24 

expressions 0(fc) below corresponds to terms of the form __I i ^ _ '^^ ''" i ^^"^^ _ '^P 

560 3 

T^pk-^ {i = 0,...,k). Consistently we have qEA{T, H) = 2/ _ 8/ _ 11/ 50/ 

^ ^27 
r + 0(2) and go(^^, 7") = P + 0(2). The computation has / 4,- 34 872 r'^' 



been carried on up to 15th order in the order parameter by 
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methods similar to those described in Ref. 12,. The only \ 27 9 / 105 

fact to take care is that having two variables to expand The expansions of energy and magnetization have been 

on, time and memory grow very fast with order, therefore computed from the following formulas [OIIMT^ 
we reduced the problem to a single-variable one consid- 



M = / P{x,y)m{x,y)dy (11) 

ering temperature-dependent values of the field, i.e. we 



= ^1 - ^ q^{x)dx^ - MH. (12) 



first computed the expansions at 15th order in r assum- 
ing j3 = r/m for m = 2, 3 ... 16, and later we reconstructed It turned out that they verify all the usual thermodynamic 
the p-dependence of the various quantities. The series are relations such as S" = —df/dT ^ and m = —df/dH. 
supposed to be valid for small r and p; although they are ^From our analysis the following conclusions can be 

non-convergent they can be resummed through the mcth- drawn regarding the PaT scaling hypothesis. 

, , . 1 . r rrm . , i . . i . , i ^ In Ref. 9 , it was sueeested that this thermodynamic rela- 

ods discussed m Ret. [12| yieldmg precise results m the _ _ 00 

tion may not be valid in the spin glass phase. The argument 

whole spin-glass phase. 

, „ „ ^ was based on an extrapolation of relatively low order results. 

^rrrr., L 2t 13t2 256 T^* \ ^ ^ 

qo{H,T) ^p IH 1 

\ 3 9 81 / It turns out from the analysis of the higher order terms that 

3/4 32 t\ 56/ ^, , , , „ 

-\-p I — — H — j — h 0(5j, (i) quantities which were assumed in [2] to depend on the temper- 

. ^ . ature alone, do have a magnetic field dependence too. These 
5t4 171 t5 1077 t6 



2 10 10 findings change the conclusions basically, leading to a support 

_ 208 t \ _ 16/ 

~^P \ 15 / 9 ^ \ ^ ' of the validity of that thermodynamic relation. 
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Fig. 1. Plot of qEA(T,H) - qEAiT,0) at fixed T = .6 as a 
function of for < < Hat{T) = .3826, the straight line is 
qsK{T, HAT{T))-qEA{T, 0) = .50689-.50504 = .00184 at T = 
.6. The data were obtained by the point-by-point procedure for 
multivariate expansions 12 using the (7, 7) Pade approximant. 
For small H, the data fit a power- law behavior .053854 ii'^"''^ 




Fig. 2. Plot of qo{H,T) at fixed = .5 as a function 
of T for < T < Tat{H) = .5351, the straight line is 
qsK{H,TAT(H)) = .5951 at H = .5. The data were obtained 
by the point-by-point procedure for multivariate expansions 
[T^ using the (7,7), (6,8) and (8,6) Pade approximants. 

— QEAiT, H) — qEA{T). This scaling law is violated in 
the expansion in power of t and p at order (i.e. 
i7io/3) which means fifth order in the height qo{T, H) 
of the first plateau. As a consequence, such a violation 
is hardly seen numerically, and the PaT approxima- 



tion qEA{T) ~ qsK{T, Hat{T)) is excellent at high 
temperature. We note that if this relation were true, 
the expansion of qEA{T) around T = should read 
QEAiT) = 1 - 3/2T2 + 0{T^). This is the expression 
which is valid on the AT line while the true expression 
is qEA{T) = 1 - aT^ + 0(7-3) ^j^j^ ^ ^ i qq ^ gi ob- 
tained resumming the expansion in power of t and nu- 
merically 112). In figuredwe plot qEA{T, H)-qEA(T, 0) 
at fixed T as a function of the magnetic field. The data 
were obtained from the series expansions by the resum- 
mation procedures discussed in Ref. '121. 

— q{x,T) — qu{j3x) is violated near the critical temper- 
ature, as can be seen from the series expansion |12| . 
However, the scaling is a very good approximation in 
a wide range of temperatures. It becomes exact in the 
limit r ^ and x/T ^ 0. A consequence of this scal- 
ing law would be \m\T-^o xi = 1/2 [S], while the real 
value turns out to be .548 ± .005 12 . From the above 
power series it can be also deduced that xi{T, H) picks 
up a small dependence on H , again proportional to 
(i.e. 

— q{x,T,H) = q{x,T). This scaling law is violated in 
the expansion in power of t and p at order p^ (i.e. 
7^10/3') -vvhich means fifth order in the height qo{T,H) 
of the first plateau. As for qEA{T), the deviation is very 
difficult to be seen numerically. Notice that this scal- 
ing law is also valid near when the more stringent 
q{x,T) — qu{x/T) is not verified. We have checked at 
15th order that all the coefficients in the expansion of 
q{x, T, H) in power of x pick up corrections propor- 
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tional to p such that the lowest power of p is always p^ 
(i.e. H^'^^^). This is likely to remain true at all orders. 
— qQ{H,T) — q{H). This scaling law is violated too. 
However, at all values of the magnetic field the be- 
havior of the function q^{H^T) is similar to that of 
figure 121 obtained by resumming the series expansion. 
Its value increases linearly from the PaT value at T = 



In particular, the field may be expressed up to fourth 
order in go hi terms of the function q{x, T, 0), we have 
indeed h"^ = 2/{3Tq'{0,T))q^ + Oqo+ 0{ql); this ex- 
plains why the scaling qo{H, T) = q{H) is best verified 
at low temperature, since Tq'{0,T) tends to the con- 
stant value 0.743 ± 0.002 for T ^ [H]. 

The predictions from the perturbative expansions have 



Tat{H) while lowering the temperature, but it goes ^^^^ compared with those from the numerical solution of 
rapidly to a constant value at low temperature. Fur- ^j^^ equations for the qix) using the high precision tech- 
thermore, at aU fields the difference between the PaT ^.^^^ introduced in ref. ^ By comparing the = 
value qoiH,TATiH)) and qo{H,0) is always of order ^-^^ the H ^ solution in aU cases we found a rather 
10-2. Actually, this is the only violation of the PaT ^^^^ agreement with analytical results and the violation 

of the PaT scaling can be appreciated. 

The series expansion at high orders allows us to safely 
sion in power of r and p which violates this additive determine the order of the transition on the AT line. By 
law is rpS i.e. rijio/^, a sixth order term. As a conse- ^^^-^^^ expansion of the rephca symmetric (RS) solu- 
quence, the deviation from the PaT entropy © is only ^-^^ ^^^^ transition line the following conclusions can 

be drawn. 



scaling which can be easily seen numerically. 
- F{T, H) = Fi{T)+F2{H). The first term in the expan 



fifth order, and again it is a very good approximation. 
The same is true for the PaT magnetization ||2Jl. We 
have checked at very high order in the expansion that 
the lowest power of p in the mixed terms is always p^ 
(i.e. i/^"/'^), i.e. fifth order in the height qo of the first 
plateau, a result that is likely to be valid at all tem- 
peratures, explaining the goodness of the PaT approx- 
imation in the whole RSB phase. Indeed, it turns out 
that the approximation q{x,T,H) = q{x,T,0) breaks 
down at fifth order in the expansion in powers of qo at 
a given finite t. Therefore we expect that at all tem- 
peratures the corrections to the PaT estimates are of 
order q^ , and become relevant only at very high fields. 



The RS-RSB transition is third order, i.e. all the pos- 
sible derivatives of the free energy (temperature, field 
and mixed) up to second order, inclusively, are contin- 
uous along the AT line. All higher order derivatives are 
discontinuous, and not divergent, at a generic point on 
the AT line. 

Near the critical point {t = Q, H ~ Q) the situation is 
more complicated. First of all, we note that the RS free 
energy is singular at this point, therefore its derivatives 
depend on how this point is approached. For instance, 
temperature derivatives of the RS free energy higher 
than the fifth are divergent approaching the critical 
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point on the AT line, while they are regular on the that at H strictly zero. Indeed at _ff = we must con- 
line {H = 0, T > 1). sider the T > 1 RS solution whose free energy differs 
Approaching the critical point (t — Q, H = Q) along from the RSB free energy at third order in r, as can 
the AT line, the discontinuity in the third field deriva- be checked using Frs = — 1/(4T) — riog2 valid for 
tive goes to zero as y^, while the fourth derivative T > 1. 

of the RS and RSB free energy (and also their dif- rpi • j i • i 

^ inese results were obtamed analyzing very long series 

ference) diverge. This fourth derivative with respect ^^^^ ^^^^^^ ^^^^^ however to confirm them we performed 

to the magnetic field is essentially the nonfinear sus- ^ ^^^-^^^ expansion of the RSB solution near the AT fine in 

ceptibility, and it can be expressed, by a fluctuation- nj. u u j-cjj- j- rpi • • 

^ jr- ' J power ot an = n ~ Hat at fixed temperature, ihis is an 

dissipation-like formula, in term of the so-called longi- • ^ li r in i u 

^ ' improvement on the result oi lU where such an expansion 

tudinal correlation function (Gl). Gl becomes iden- • j • men -c j-- • ^ j -cj-i. i 

^ ^' ^ was obtained using a IRbB function instead of the correct 

tical with the spin glass susceptibility in exactly zero r ,1 octd ^ \ rru- • n 4- j- j-i 

^ ° J J j-^ll KSB q[x). ihis expansion allows to compute exactly 

magnetic field, and diverges at T, (r = 0), although it derivatives of the RSB free energy on the AT fine, 

remains finite elsewhere. expected we have found that the first and second field 

The discontinuities in the third and fourth tempera- j ■ i- r 4-\ -do j octd t ■ i- 

^ derivatives of the RS and RSB free energies are contmu- 

ture derivatives also go to zero while approaching the iUATni- t 1 u j ^i. ^ ■ 

° e, ous on the AT line. In general we have observed that m 

critical point along the AT line, only the fifth temper- j^ji-j-if jun ■ rc 

f o t J f order to determine the tree energy at order an it is surfi- 

ature derivative has a finite discontinuity in this fimit. ^.^^^ ^^^^ ^^^^ ^^^^^^^ ^^„„2 -^^^^^^ ^f order dh"'^ 

As already mentioned, higher temperature derivatives ^^-^^ ^^^j^ expected by power counting. Therefore 

of the RS free energy diverge at Tc, while the RSB ^j^^ solution, which is the zeroth order approximation, 

derivatives are regular. Therefore, if we consider only • .1 . r , j j • jj, 

' gives the correct tree energy up to second order m an,, ihe 

temperature derivatives, the RS-RSB transition be- eoefiicients of this expansions depend on integrals of the 

comes fifth order near the critical point. This can be r 

^ form 

understood noticing that, while approaching the criti- -, ^-1-00 

^=-=/ e-y'^'^^'^ha.iih[py + pHYdy (13) 
cal point from the left along the AT line, the RS solu- V^^? J -00 



tion, valid above the line, is similar to the (unstable) 



evaluated on the AT line. For instance the second field 



RS solution at T < 1 and = whose free energy derivative reads 
differs from the RSB free energy at fifth order in t. — — — ( T Hat(T)) 



dH^' ' " 1 + (-7AT - l)/52 

This behaviour of the temperature derivatives on the / , r> r)2 2 ^2 

X [qAT - 1 - 2qAT/3 + qArP 

AT line in the iJ — > limit is markedly different from 2 



-0^{-l + Mlj,~2MATm^ + mf)^ (14) 



8 
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The breakpoint reads 



13. H.-J. Sommers, W. Dupont, J. Phys. C 17, (1984) 5785. 



2(-2 + 2/32-3gAT/3' + /3'm6) 
3 - 2/32 + SqATfi^ - P'h^ 



14. A. Billoire and B. Coluzzi, Phys. Rev. E 67, (2003) 036108. 



xx{T,Hat{T)) 



(15) 



It is interesting to notice that the T ^ limit of the 
breakpoint is exactly 1/2 as predicted by the PaT hy- 
potheses in an independent way and at variance with the 
(T = 0, = 0) value which is .548 ± .005 [H]. 
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